Angular position of nodes in the superconducting gap of YBCO 
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Abstract 

The thermal conductivity of a YBa2Cu30g.9 detwinned single crystal has 

been studied as a function of the relative orientation of the crystal axes and a 

magnetic field rotating in the Cu-0 planes. Measurements were carried out at 

several different temperatures below T c (0.5 K < T< 25 K) and at a fixed field 

of 30 kOe. A four-fold symmetry characteristic of a superconducting gap with 

nodes at odd multiples of 45 degrees in k-space was resolved. Experiments 

were performed to exclude a possible macroscopic origin for such a four-fold 

symmetry such as sample shape or anisotropic pinning. Our results impose an 

upper limit of 10% on the weight of the s-wave component of the essentially 

d-wave superconducting order parameter of YBCO. 
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Over the last years, there has been accumulating evidence that the pairing state of 
high temperature superconductors may be the so-called d x 2_ y 2 state. Among the wealth 
of experimental data, one can mention the results from magnetic penetration depth 
angular resolved photoemission and nuclear spin relaxation rate |J, which provide strong 
evidence for nodes in the superconducting gap of the cuprates. Furthermore, quantum phase 
interference experiments |§ which look directly at the symmetry of the order parameter as 
a function of its argument, generally confirm the presence of the nodes and the sign change 
of the gap function over the Fermi surface. 

Most of these experiments suggest, as a recent review by Annett et al. f7| concludes, that 
a d x 2_ y 2 order parameter is the most plausible candidate to describe the superconducting 
state in these systems. However, in the case of optimally-dopped YBa2Cu307_,5 (YBCO), it 
has been argued that a possible slight s-wave component could arise from the orthorhombic 
distortion ||. The in-plane anisotropy observed below T c in thermal conductivity and 
penetration depth suggest that the CuO chains play a role in the thermodynamic and 
transport properties of the superconducting state. Furthermore, it was recently pointed out 
that the in-plane penetration depth anisotropy is difficult to understand within a proximity 
model with only a pure <i-wave pairing in the Cu-0 plane |ItJ. One consequence of the 



existence of a s-wave component would be a shift of the nodal positions from the principal 
diagonal directions over the Fermi surface. Thus, the precise angular location of the gap 
nodes on the Fermi surface would be a useful information in the current debate on the 
symmetry of the order parameter in YBCO. 

We present here a study of the thermal conductivity of a YBCO single- crystal, k, ro- 
tating in a magnetic field parallel to the Cu-0 planes. A pioneer work by Salamon and 



co-workers [fTT] , [T2|| showed that thanks to the Andreev scattering of quasiparticles by ro- 
tating vortices, such an experiment can be a relevant probe of the k-space anisotropy of 
the superconducting gap. This work was criticized by Klemm et al. JTB) who contested a 
microscopic interpretation of the data. Moreover, error bars on the angular dependence of 
the longitudinal thermal conductivity were comparable to the magnitude of the periodical 



oscillations [pj] , making it difficult, if not impossible, to resolve a four-fold symmetry of k{6) 



which is the qualitative distinction of d-wave superconductivity. In this letter, we begin by 
reporting on high-resolution experimental results which establish such a four-fold symmetry 
and then we address objections to a microscopic interpretation of this angular dependence. 
Finally, we show that our data impose an upper limit on the relative weight of the s-wave 
component of the essentially d-wave state of YBCO. 



Our YBCO single crystal was prepared by a self-flux method described elsewhere |14 
Measurements were performed in a dilution refrigerator using a one heater-two-thermometer 
steady state method. The sample was anchored at one end to a Cu heat sink, and a heater 
resistor was attached to the other. The heat current (J) was injected along the b-axis 
direction (i.e. parallel to the chains), and the thermal gradient was measured along the same 
direction with two Lakeshore Cernox thermometers. The experimental setup was rotated 
at the center of a superconducting coil with the rotation axis along the c-axis and the 
magnetic field parallel to the Cu-0 planes. The precise relative orientation of the magnetic 
field and the crystal was determined through a Hall probe rigidly fixed on the setup. The 
maximal available rotation was about 200 degrees. Below 5 K, due to the slightly anisotropic 
magnetoresistance of the thermometers, we used a dynamic method which allowed us to 
calibrate the thermometers for each measured ^-direction. 

Fig. 1 shows the thermal conductivity of the YBCO sample as a function of temperature 
over a wide temperature range. A few years ago, it was suggested that the steep increase 
of the thermal conductivity at the superconducting transition is induced by the electronic 
contribution in the Cu-0 planes which, due to the strongly suppressed quasiparticle scat- 
tering rate in the superconducting state, increases rapidly below T c 0. Early evidence for 
such a suppression came from microwave conductivity data Later, thermal Hall effect 
measurements (16| confirmed this pattern and it is now generally accepted |I7j that the peak 
of the thermal conductivity below T c is in large part due to electrons. In our sample this 
peak is higher {^f^ = 2.2) and occurs at a lower temperature (~ 23 K) than what has been 
previously reported [[J. Both these features point to a long maximum quasi-particle mean- 



free path in the sample studied in this work. Because of this large electronic contribution, 
thermal conductivity is a useful probe of the quasi-particle excitation spectrum. 

The angular variation of the thermal conductivity of YBCO, k (6) with 9 = (b, H), at 
24 K is shown in fig. 2. The full range of angular variation ( < 6< 360 degrees ) was 
obtained by inverting the direction of the magnetic field. Thermal conductivity shows a 
clear fourfold variation superposed on a hysteretic background. This hysteresis, observed 
between successive rotations, is related to the pinning of vortices, which might induce some 
variation of the vortex density between two consecutive crossings of the magnetic field at a 
^-direction. This hysteresis continues to be present at lower temperatures where the pinning 
is stronger. But, somewhat surprisingly, the monotonic background disappears. As shown 
in fig. 3a, the angular variation has a clear fourfold symmetry at 6.8 K. This is the first time 
that such a symmetry of thermal conductivity is unambigously established well beyond the 
experimental resolution in a high-T c superconductor. 

Next, we consider possible objections to a microscopic interpretation of this fourfold 
symmetry such as raised by Klemm et al. [|T3]. These authors suggested that the angular 
structure might be "due to the demagnetization and flux pinning effects associated with a 
rectangular sample" . To check this hypothesis, we performed complementary verifications. 
First, we used an alternative experimental procedure by rotating the crystal in the normal 
phase: After each rotation of about 1 degree, the sample was heated to a temperature above 
T c (93 K) before cooling down to the measurement temperature. The angular variation 
obtained in this way ( fig. 2) presented the same behavior without the sloping background. 
Second, to rule out the possibility of a morphological origin -related to the demagnetization 
for example-, we studied carefully the thermal conductivity of a niobium (Nb) sample of 
similar shape (flat rectangular slab) and size and our results reproduce previous measure- 
ments on Nb [[TjJ. As shown in fig. 3b, the angular variation of the thermal conductivity for 
this sample has an obvious twofold symmetry. This anisotropy changes sign with decreas- 
ing temperature reflecting the nature of the dominant quasi-particle carriers [jn5|, but the 



symmetry always remains twofold. Thus, the relative orientation of vortices and the applied 



current, and not the sample geometry, governs the angular structure seen in fig. 3b. 

Now, with the elimination of a macroscopic origin for the fourfold variation observed in 
the case of YBCO, we can turn to the microscopic picture which was originally formulated 



by Yu et al. (12| and is based on the anisotropic scattering of the quasi-particles by the 
vortices. It is well-known |2D| that the vortices have dramatic effects on the heat transport 
in type II superconductors. Even in the case of the cuprates with a magnetic field parallel 
to the basal plane, the coreless vortices lying between the Cu-0 planes can act as scatterers 
of the BCS-like quasiparticles. This scattering mechanism was expressed || in terms of 
Andreev reflection |HJ] which is a process of retroreflection of excitations where spatial 
variations of the amplitude or the phase of the order parameter induce branch conversion of 
electron-like excitations into hole-like excitations, and vice versa. In other words, the excited 
states distributed above the superconducting gap exhibit in the presence of a phase gradient 
(superfluid flow, v s ), a Doppler shift of their energies, E=E — p.v s , so, when a quasiparticle 
approaches the vortex core, its energy in the superfluid frame reaches the energy gap and 
is converted to a quasihole, halting its contribution to the heat transport. This scattering 
process is strongly dependent on the relative orientation of the quasiparticle momentum (p) 
and the superfluid velocity (v s ). The latter is imposed by the magnetic field. Indeed, no 
Andreev reflection occurs when the quasiparticle momentum is normal to the supercurrent 
and parallel to the magnetic field. Thus, thanks to the directionality of this scattering 
mechanism, thermal conductivity can identify the preferential momentum orientation of the 
excitations above the superconducting gap. 

Qualitatively, the observed behavior, i.e. a fourfold variation with maxima for odd multi- 
ples of 45°, indicates that there are maxima in the angular distribution of the quasi-particle 
momentum in the vicinity of | k x | =| k y | . When the field is aligned along these particular 
directions, most quasiparticles are not Andreev reflected, and the thermal resistance is lower. 

The predominance of the fourfold variation at 6.8 K enables us to compare the curve to 
a 2D version of the usual BRT expression of the electronic thermal conductivity |T3|: 
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where v g b is the the b-axis component of the group velocity, c is the c-axis lattice parameter 
and the effective rate of scattering 



| ,tt 1 M 4- T ° Pvn I ~ m2Q " [^P- I A (P)U I , 
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depends on a v , the intervortex mean distance, and on i/>(p), the angle between the mag- 
netic field and the quasi-momentum direction p. — and — are respectively zero-field and 
maximum Andreev scattering rates. The gap function is a mixture of d-wave and s-wave 
components: 

A-A [r I cos ^ a /^) ~ cos (Py a /h) 1 
\ 1 — cos{pFd/Ti) J 

where r represents the relative weight of the s component. As mentioned above, a finite r 
would shift the nodal positions from the main diagonal directions and the angular interval 
between the peak positions of the thermal conductivity would differ from 90°. Numerical cal- 
culations using realistic physical parameters (E F = leV, r = 3r VQ = 10 _12 s and — = 0.05) 
are shown in fig. 3a for pure d — wave (r = 0), d + 10%s(r = 0.1) and d + 30%s(r = 0.3) 
cases. The interval (29) between the two peaks is shifted from 90 degrees (pure d — wave) 
to 96 (10 %) and 108 (30%), following the relation 28 = Arccos(r). However, in our experi- 
mental curve, the two peaks are separated by 90±6°. Thereafter, according to our results, 
the upper limit to r, i.e. the relative weight of the s component of the superconducting 
order parameter, (A = A^(r + cos(2#))), is about 0.1. This is significantly lower than what 
is suggested through a recent examination [^2| of the in-plane penetration depth anisotropy 
0]. On the other hand, it is compatible with the absence of any detectable difference in 
the zero-field thermal conductivity along the a-axis and the b-axis at very low temperatures 
which also suggest an essentially d-wave state |23| . 



We also studied the temperature dependence of this angular variation down to 0.5K. 
Fig. 4 shows the curves for 24K, 6.8K and 0.8K plotted together. A striking change in the 



structure of k{9) is detected between 6.8K and 1.2 K where the fourfold symmetry fades 
away. The minimum for the field perpendicular to the thermal current (b-axis) develops 
significantly (i.e. Kn > k±). We have found that this minimum becomes deeper with 
decreasing temperature or for increasing magnetic field. The origin of this crossover between 
these two regimes (6.8K and 1.2 K) is not yet well understood. 

However, the analysis above neglects the effect of vortices in the plane on both phonons 
and quasiparticles in the chains. The former may provide an explanation for the twofold 
variation developed at low temperatures. On the other hand, as a careful study of the in- 
plane anisotropy of thermal conductivity at various temperatures suggests ||24|| , the chain 
contribution to the thermal conductivity increases significantly below 50K passing by a 
maximum around 15 K. In our geometry, with the heat current applied along the chain 
orientation, a two-fold angular structure due to the scattering of chain quasi-particles by 
vortices could well be present. Additional experiments of this type on a-oriented YBCO 
crystals will be necessary to elucidate the effect of a rotating magnetic field on different 
types of heat carriers at different temperatures. 

In conclusion, our results support the existence of nodes in the superconducting gap at 
angular positions close to what is expected for a purely d-wave state and impose an upper 
limit of about 0.1 to the relative weight of the s component. 
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FIGURES 

FIG. 1. Main Panel: The temperature depencence of the thermal conductivity in a log/log 
plot. Insert shows a linear plot. Note the sharp increase at T c leading to a peak at 23 K. 

FIG. 2. Open squares and circles show the angular variation of the thermal conductivity at 
T=24 K and H= 30 kOe. There are minima for the field parallel and perpendicular to the b-axis 
direction (heat current), and maxima when the field is parallel to the main diagonal directions. 
Arrows indicate the rotating direction. Solid circles indicate the results obtained with an alternative 
procedure of field-cooling at every angle (see text). 

FIG. 3. The fourfold variation of the thermal conductivity in YBCO (a) is compared with the 
twofold variation for the Nb crystal (b). For YBCO the numerical results are shown for a pure 
d-wave, d + 10%s — wave(r = 0.1) (solid line) and d + 30%s — wave(r = 0.3) (dashed line). For 
Nb the solid line is a fit using the cosine-square dependence expected for a conventional gap. 

FIG. 4. Angular variation of the thermal conductivity at 24 K, 6.8 K and 0.8 K, at H= 30 kOe. 
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